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ABSTRACT 


A general  method  is  presented  for  obtaining  the  dynamic  equations 
of  an  elastic  structure  to  which  elastic  and/or  nonlinear  substructural 
elements  are  attached. 
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I INTRODUCTION 

Component  mode  methods  had  been  successfully  exploited  to  solve 
vibration  problems  long  before  the  computer  age  made  it  feasible  to 
analyze  complex  structural  systems  by  discretization,  e.g.,  the  finite 
element  method. 

While  analytical  solutions  abound  in  the  literature,  the  structures 
considered,  out  of  necessity,  are  comparatively  simple,  and  are 
analyzed  on  an  ad  hoc  basis.  The  extensive  use  of  finite  element  methods 
has  made  it  feasible  to  consider  the  most  complex  structural  systems, 
and  thereby  made  it  necessary  to  develop  comprehensive  procedures  for 
obtaining  the  dynamic  equations  for  the  system. 

Various  methods  for  component  mode  analysis  have  been  suggested  [1,  2,  3]. 
They  differ  from  one  another  in  their  choice  of  component  modes,  but  all 
express  the  equations  of  motion  in  terms  of  system  modal  stiffness  and 
mass  matrices  and  component  mode  coordinates.  Their  concern  is  essentially 
directed  toward  solving  the  eigenvalue  problem  and  determining  the  system 
modes . 

In  a previous  report  [4]  a general  method  was  presented  for  obtaining 
the  dynamic  equations  of  an  elastic  structure  to  which  arbitrary  elastic 
substructural  elements  are  attached.  The  methodology  used  was  to  express  the 
kinetic  and  potential  energies  of  the  system  as  the  sum  of  the  contri- 
butions from  the  main  structure  and  substructures,  separately.  It  was 
shown  that  the  equations  of  motion  of  the  main  structure  were  affected  by 
the  presence  of  the  attached  structure  by  means  of  the  forces  exerted 
at  points  of  attachment. 
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It  is  the  purpose  of  this  report  to  formalize  and  fully  automate 
the  procedure,  and  to  extend  it  so  that  it  can  be  applied  to  the  most 
complex  structural  elastic  systems  having  nonlinear  attachments  and/or 
small  substructural  regions  that  respond  nonlinearly.  A treatment  of 
some  special  cases  involving  nonlinear  attachments  may  be  found  in  [5]. 
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11  GENERAL  METHOD 

4 

Consider  an  elastic  structure  S to  which  an  elastic  substructure  O is 

attached.  We  denote  the  in  vacuo  modal  and  diagonal  generalized  mass  matrices  of 

S by  <t>g>  )Jg,  respectively,  and  the  corresponding  diagonal  frequency  matrix 
2 

by  Wg,  and  express  the  displacement  of  the  structure  as 

dS  = Vs  (1) 

where  qg  is  the  generalized  coordinate  vector. 

The  kinetic  and  potential  energies  of  the  system,  expressed  as 
sums  of  contributions  from  S and  0 separately,  are 

T(V  V = 2 % Vs  + 2 Wo  ^ 

v (qS , da)  = | qg  ysWgqs  + \ d^K^  (3) 

where  and  are  the  mass  and  stiffness  matrices  of  the  substructure, 
t and  dQ  is  the  0-displacement  vector. 

c 

Let  the  constrained  base  modes  be  defined  as  the  static  elastic 
and  rigid  body  displacements  of  0 due  to  successive  unit  displacements 
of  the  constrained  degrees-of-freedom  (dof)  (i.e.,  those  interface 
coordinates  of  0 which  are  constrained  to  move  with  S).  Define  the  fixed- 
base  modes  as  the  modes  of  0 obtained  when  the  constrained  dof  are 
held  fixed. 

The  displacement  of  the  substructure  is  expressed  as  a superposition 
of  its  static  response  due  to  movement  of«  its  constrained  degrees  of 
freedom  and  its  motion  relative  to  them; 
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This  can  be  expressed  in  a most  convenient  compact  form  [3]  as 


1 , 
.^D  ! 
I l 

= 

<j>  1 

rS  ' 

1 

0 

_3°J 

j 

TC6  1 

L s i 

A 

^fo_ 

*fa 


(5) 


where  ( ) and  ( ) refer,  respectively,  to  the  constrained  and  unconstrained 


dof;  T is  the  constrained  base  mode  matrix  whose  successive  columns 


are  the  constrained  base  mode  vectors;  q represents  the  fixed-base 


generalized  coordinate  vector. 

For  convenience  we  express  the  displacement  of  Eq . (5)  as: 


c*— 

d = T d q„  + <£„  q 
a yS4S  vfo4fo 


(6) 


where 


(7) 


and 


Kf0 


1-0- 

$ 


fa 


(8) 


Lagrange's  equations  for  S are 


_9V_ 

dt  3V 


(9) 


where  Qs  are  the  generalized  forces  due  to  external  forces. 


Substituting 


Eqs.  (2)  and  (3)  into  (9)  then  yields 


2 3dT 

Wa  + Ws  + Vo  + Koda  = Qs 

b b 


(10) 


But  from  Eqs.  (6) 


3do  “a 

aq  ' m;  ’ *s  T 


(ID 
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and  thus  Eq.  (10)  can  be  written  as 


2 -T  c*T 

M q + M a)  q + <J>  T (M  d + K d ) = Q 

J b bb  b b O Q UU  b 


Similarly  by  applying  Lagrange's  equations  to  the  substructure  we  obtain 


M d + K d = F 

a o a a 


If  we  make  the  assumption  that  no  external  forces  are  applied  to  the 
substructure,  the  force  vector  F vanishes  everywhere  except  at  the 
constrained  dof.  Thus 


where  F are  the  interaction  forces  exerted  on  a at  the  constraints. 

Note  that  no  difficultes  arise  when  external  forces  applied  to  the  sub- 
structure are  included  (see  Appendix  A) . 

Upon  substituting  Eqs.  (7),  (13)  and  (14)  into  (12),  we  can  write 
the  equations  of  motion  for  the  main  structure  as 

2 -T  - 

p q + p u)  q + <{>  F = Q 
b b b b S S S 

The  equations  for  the  substructure  (13)  can,  by  substituting  Eqs.  (6-  8, 
14),  be  recast  as 


“oVs  + % + Vs  + Vfaqfa  = 5 


(12) 


(13) 


(14) 


(15) 


(16) 


Vfaqfo+  Vfo^fa  + VC  Vs  + (*a  + Vs  = 0 (17) 


where  m,  k are  the  elements  of  the  partitioned  mass  (assuming  a lumped  mass 
approximation)  and  stiffness  matrices  : 
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m l 

o 1 

fk 

k 1 

a | 

; k = 

’ a 

_q_  o 

A 1 

0 

m 

£ 

k 

L ‘ aJ 

a 

aJ 

(18) 


By  the  very  definition  of  the  constrained  base  modes  [i]  it  follows  that 


<ka  + k0T  ) <t>sqs  5 0 


(19) 


Consequently,  upon  pre-mult iplying  Eq.  (17)  by  4>f0  , we  have 


tJfoqfa+  wfo“faqfo  + hi  ("aT%5s)  = 0 


(20) 


where  p and  represent,  respectively,  the  generalized  mass  and 

frequency  matrices  of  cr. 


In  summary,  the  dynamics  of  the  system  is  expressed  by  Eqs.  (15) 
and  (20), where  F,  the  interaction  force  vector  is  given  by  Eq.  (16). 

Note  that  the  first  term  of  Eq . (16)  is  that  due  to  the  inertia  of  the 
mass  of  0 located  at  the  constrai ned  dof , while  the  remainder  of  the 
expression  is  due  to  the  stiffness  forces  exerted  by  a (see  figure  below)* 

' > 

lnlp-1 

Tf 

S_  

INTERACTION  FORCE 


There  are  a number  of  advantages  in  considering  the  main  structure 
separately,  and  accounting  for  the  substructure  by  means  of  the  constraint 
forces.  First,  the  computer  core  (central  memory)  required  to  solve  a 
given  problem  is  reduced,  since  no  total  system  modal  stiffness  matrix  is 
needed.  Second,  once  the  constraint  modes  have  been  determined,  only 
those  rows  of  (the  unconstrained  stiffness  matrix  of  o)  that  correspond 


7 


| 


1 

i 


V 


to  the  constrained  physical  degrees  of  freedom  of  o are  required  in  the 
solution.  This  also  reduces  memory  requirements.  Third,  the  amount  of 
computer  time  (central  processor  time)  required  to  set  up  the  matrices 
in  the  equation  of  motion  is  reduced,  as  compared  to  that  required  by 
a method  employing  a system  modal  stiffness  matrix.  Referring  to  the  method 
of  Ref.  [3],  for  example,  the  system  modal  stiffness  matrix  is  determined 
by  the  second  of  Eqs.  (A-2) . The  computer  time  needed  to  evaluate  this 
matrix  may  be  excessive,  especially  if  the  problem  under  consideration 
requires  many  modes  and/or  physical  degrees  of  freedom,  and  if  is  not 
well  banded.  In  fact,  the  rearrangement  of  physical  degrees  of  freedom 
of  0 to  partition  usually  results  in  a poorly  banded  matrix.  Similar 
difficulties  with  the  first  of  Eqs.  (A-2)  may  be  avoided  by  using  a 
lumped  mass  approximation  for  o. 

In  view  of  the  above,  the  method  developed  here  allows  one  to  more 
readily  handle  complex  systems  with  many  substructures  within  the  re- 
strictions imposed  by  a given  size  of  central  memory. 


FREE-FREE  MODES  OF  SUBSTRUCTURE 


The  displacement  of  the  substructure,  expressed  in  terms  of  its 
free-free  modes  $ , is: 


d 

0 

l 

<j>„  i 0 
s ! 

' \ ' 

— 

r" 

d 

0 : * 

q 

0 

i 0 

_ 0 _j 

(21) 


since 


*0  ^ = h*S 

a 


Proceeding  as  before,  we  obtain  the  following  set  of  equations  in  lieu 
of  Eqs.  (16)  and  (17). 


(22) 


r - 


•>  . 
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I 


1 

\ 

i 


m <p  q 

a oua 


+ Volo 


k 4>  q 
oYo^a 


= F 


(23) 


ma*oqa 


ko*c% 


ka4’oqa 


= 0 


(24) 


The  generalized  mass  (p^)  and  frequency  (oj^)  matrices  for  the  free-free 
modes  of  the  substructure  are 


X — T—  — ^ 

u = d)  M cb  = d>  m (b  + d)  m 4> 

Ma  ya  oYo  ya  era  Ya  oyo 


u cj2  = (f^K  d>  = d)"^  (k  4>  + k d>  ) + <l)^(k  <J)  + k (|)  ) 

Mo  a Ta  oyo  Ta  a a aya  a era  oyo' 


(25) 


If  we  now  pre-multiply  Eq.  (24)  by  <{>  , and  make  use  of  the  expressions 

in  Eq . (25),  then  the  following  is  obtained 


Vo  + Wo  ' + (Vsqs  + kaVa)1=  0 


(26) 


The  motion  of  the  main  structure  is  determined  by  Eq.  (15),  which 
when  combined  with  Eq . (23)  becomes 


Vs  + Ws  + (Vsqs  + Vaqa)]=  Qs  (27) 

Thus,  when  free-free  modes  of  the  substructure  are  used,  Eqs.  (26) 
and  (27)  express  the  motion  of  the  system.  Note  that  the  bracketed 
terms  in  (26)  and  (27)  represent  the  constraint  interaction  forces, 
and  that  the  motion  of  the  substructure,  Eq . (26),  is  coupled  to  that  of 
the  main  structure  through  the  enforcement  of  the  constrained  dof  relationship, 
Eq.  (22). 

NONLINEAR  ATTACHMENTS 

Consider  a composite  structural  system  in  which  each  of  the  components 
is  assumed  to  exhibit  a linear  elastic  response,  and  in  which  the  components 
are  attached  by  nonlinear  mountings.  We  denote  the  force-displacement 
relationships  of  the  mountings  by  F(6,  6),  where  6 is  the  relative  inter- 


face displacement  vector 


/ 


— 9 — 


6 = ds  ' da 


(28) 


When  component  mode  analysis  is  to  be  applied,  it  is  necessary  to 
isolate  these  attachments  from  all  of  the  components,  so  that  only  the 
forces  transmitted  by  them  enter  into  the  analysis. 

FIXED-BASE  MODES.  The  displacements  of  0 are  expressed  as 


a l n 

s _J 

TCd)  6 

^S  I vfa 


1 


L£0j 


(29) 


When  determining  the  fixed-base  modes  of  O,  it  should  be  noted  that  the 
constrained  dof  must  include  those  which  are  constrained  by  the  nonlinear 
attachments. 

Upon  combining  Eqs.  (13),  (14)  and  (29),  we  obtain  the  following 
equations  for  the  substructure1 


and 


F = m ($_qc  - 6)  + (k  + k TC)  (<j>cqc  - 6)  + k 4>_  q. 


Ufoqfa+  wfo“fa  qfo+  ♦fo  VC(*SqS  “ 6)  = ° 


(30) 


(31) 


The  nonlinear  constraint  forces  are 

F = F(6,  6) 

Thus  the  system  dynamics  is  expressed  by  Eqs.  (30),  (31),  (32)  and  (15). 

FREE-FREE  MODES.  When  free-free  component  modes  are  used,  we 
express  the  displacements  as 

I /s 

,°  1 +a 


(32) 


(33) 
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Utilizing  Eq.  (33)  and  proceeding  as  before,  we. obtain 


F = mo($sqs  - 6)  + k(iq5  - 6)  + k Ja 


yoqo  + yoVo  = *oI“o(*siis  " g)  + V^s  ’ 6)  + Vaqo] 

The  substitution  of  Eq.  (34)  into  (35)  results  in 

%%  + Wo  = *o  «> 

where 


6 = Vs  " Vo 


(34) 

(35) 


(36) 


(37) 


Equations  (15),  (32),  (34)  and  (36)  characterize  the  dynamics  of  the 
system  when  free-free  inodes  are  used. 

One  might,  upon  juxtaposing  Eqs.  (15),  (30)  and  (31)  and  Eqs.  (15), 
(34)  and  (36),  arrive  at  the  conclusion  that,  where  nonlinear  attachments 
are  to  be  considered,  it  is  more  convenient  to  use  the  free-free  modes 
of  the  substructure.  Certainly,  their  use  eliminates  the  need  to  store 
any  of  the  mass  and  stiffness  elements  of  component  0,  as  would  be  the 
case  when  fixed-base  modes  are  used.  Also,  only  two  coupled  systems  of 
equations  result  when  they  are  used.  However,  these  advantages  can 
readily  be  offset  by  the  need  to  use  a far  greater  number  of  free-free 
modes  than  would  be  required  when  fixed-base  modes  are  used.  For  example, 
for  rather  stiff  mountings,  fewer  fixed-base  modes  would  be  required. 
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III  ILLUSTRATIVE  PROBLEMS 


A.  Two-Degree-of-FreedOm  System  - Linear  Springs 


( i)  Fixed-Base  Modes 


k -k  0 0 

K = Z ; H = 

0 , , o „ 

-k2  k^  0 


Hfs_lis  = d<y=qs 


and  thus 


k0  " k2;  ko  = ko  ' -k2 ’ ko  “ k2 


:©  <ps  = I 


ma  = 0;  rao  = m2 


TC  = l;  $fa  = 1;  Ufa  = m2;  Ufa  = k2/m2 


Upon  substituting  into  Eq.  (20),  we  obtain  the 


1 df  cr  = ( 


*2  4>  f«r  = ' 


equation  of  motion  for  o: 


FIXED  -BASE  MODEL 


V-f,;  + k2<lfa  + m2^  = 0 (1) 


But  from  Eq.  (5):  d(j  = qg  + qfQ,  and  tnut 


-2d0  + k2<du  - 3S>  ' 0 (I  > 


For  the  main  structure 


Ms  = rai;  ws  = Vv  h = 1 


and  thus  Eqs.  (15)  and  (16)  yield: 


V-S  + Vs  ~ k2qfO  “ fs 


mlds  + klds  ~ k2 ^dO  ~ V “ ^ g ^ 
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(ii)  Free-Free  Component  Modes 


ko  “ k2;  ka  ko  "k2;  ka  " k2; 


mo  = 0;  mo  = m2 ’ 


= 


1 1 

1 0 


m2  0 
0 0 


; “ol  = °;  ">02  = °° 


Vo  = *Wo 


0 0 

0 k„ 


Us  = m1 ; 


From  Eqs.  (21)  and  (22), 


A ✓v 


do  = dS  = Vo  = Vs;  do  = Vo 


? 


\ 


and  thus, 

do  = qol  + qo2;  do  = qol 
Consequently,  Eqs.  (26)  and  (27)  yield: 

(iii) 
(iv) 


m2d0  ~ k2(do  " V = ° 


mld0  + kld0  + Vdo  ■ V = fS 


*) 

When  determining  the  free-free  modes  of  the  O-structure,  one  might, 
at  first  glance,  conclude  that  the  constrained  dof  need  not  be 
considered,  since  no  mass  is  placed  there.  But  the  o-structure  has 
two  degrees  of  freedom  and  thus  two  free-free  component  modes  are 
required.  The  first  of  these  is,  of  course,  the  rigid  body  mode 
which  has  a zero  frequency;  the  second,  which  is  not  so  apparent, 
is  a mode  in  which  the  constrained  dof  oscillates  with  an  infinite 
frequency,  so  that  the  unconstrained  dof  remains  immobile.  Alternately, 
one  can,  by  placing  a fictitious  mass  m at  the  constrained  dof, 
determine  the  free-free  modes  in  the  limit  as  m + 0. 


- r 'j. 
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(ii)  Free-Free  Component  Modes 


Ko  = 0;  mo  = Mo  = m2:  “a  = 0; 


MS  = ml;  WS  = kl/ml:  = ^ *0  = 0 


Equation  (22)  yields:  d^  - 

and  from  Eqs.  (15)  and  (36),  we  have 


fd<r 

NONLINEAR 
(£)  " SPRING 

Vs JdS 


mids  + kids  + F = fs 

(iv) 

m2da  * F = 0 

(v) 

where  F = F(6,  6)  and  6 = dg  - dfJ 


C.  Simply  Supported  Beam  and  Spring  Mass  System  - Linear  Springs 
(i)  Fixed-Base  Modes 

r~  — » I — l 1 L C 


kx  -ki 


Ko  ' 'kl  <ki+k2)  "k2  • Ka  ■ 0 »!  0 


0 -k„ 


0 0 0 


0 0 m, 


© ♦«, 


Let:  k^  = 2k;  - k;  m * 2m,  m2  = m 

Then:  ka  = 2k;  ka  = £-2  o|  k 

K {o  ] k‘  k„  ■ [-!  !]  k 

"o  ' 0;  "a  ' [o  l]  ■ 


M. 


rSl  g-fr 

Kuf 

BEAM  CONSTANTS 


S"1"'  ti'-  A = cross-sectional  area 

•Jdfai 

El  = bending  rigidity 
► . p = mass  density 


L J FIXED- BASE  MODEL 

It  can  readily  be  shown  that  the  fixed-base  modal  and  frequency  matrices 
are:  r r . -i 


1 1 
2 -1 


1/2  0 
0 2 


--  15 


and  that  the  constrained  base  mode  matrix  is: 

-Cl 

rv,  r 

I T 


Consequently,  uf(J  = <J>f*  ma^fa= 


6 0 
0 3 


For  the  main  structure: 


4 4 

pA£  2 it  n El 


, i-  n n ti  , v . liiTx 

Ps  = T;  “s  = 74-  • Tk  ' = l %n  sln  — 


Substitution  into  Eq . (20)  yields: 


,1k  . 2 c ••  . nTT  n 

*fal  + 2 m qfal  + 3 1 qSn  Sln  l = 0 


*fa2  + 2 m qfo2  + 3 ^ qSn  sin  = ° 


Now  from  Eq . (16):  F = -2k(qfQl  + qfa2> 


and  from  Eq . (15) : 


4 4 

. ti  n El  .2  .nr 

1,,-  + — 1~  • 77  + — sm 


*■80  ^4  ‘ pA  -Sn  pA£ 


-2k(afai  + q.fa2) 


= Q, 


Sn 


The  following  identities  are  established  from  Eq . (5): 


^ol  = ^fol  + ^02  ; 2 ■ 2q-fol  “ 


fo2 


Thus  (i)  + (ii)  and  2(i)  - (ii)  yield 


2mdal  + 2k(dol  - d0)  + k(d01  - da2)  = 0 


md  + k(d  „ - d , ) =0 

a 2 o2  al 


(i) 


(ii) 


(iii) 


(i’i 

(ii’) 


-£ 
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while  (iii)  can  be  written  as 


4 4 

^ it  n El  ^ 2 . Tin 

^Sn  ^4  pA  ^Sn  + pA?.  S1°  2 


where  d = d = /,  q.cnsin 


. HIT 


-2k(doi  - d0)  = QSn 


(iii') 


(ii)  Free-Free  Component  Modes 


The  free-free  component  modal  and  frequency 


matrices  are: 


1 1 1 

110 
1 -2  0 


3 0 0 


J “a  = 


0 6 0 ; p a) 

ooo  a ° 


0 0 0 

0 3/2  0 k/m  ; 
0 0 °° 


0 0 0 “ 
0 9 0 

0 0 2 


From  Eq . (26) 


V - 2 mf  l qSn  sin  f - (qol  + W]=  0 

^2  m ^2  ‘ 3 ^Sn  sin  f ~ <%1  + W]=° 


Equation  (27)  yields: 

4 4 

. v n El  . 2 . nu  r r . irrr  , , ... 

V +—  * qSn  + ^ sln  T 1 ^ V1"  T " %1  %2)]  = QSn  (vl) 


But  from  Eq.  (21) 


%1  + %2  5 do2  ‘krl  ?<*o2 


and  thus  (iv)  + (v)  and  (iv)  - 2(v)  yield 


2"  3.1  + «dol  - So2>  •2kp’sn*1"  f - 5ol]  ■ 0 «“') 


“*02  • k<d01  ‘ do2>  ' 0 


17  — 


while  (vi)  becomes 


s +»V 


El  2 nff  r c nir  P ■, 

PA  pa£  Sln  T [^qSn  sin  T " dal]  = 


D.  Simply  Supported  Beam  and  Spring  Mass  System  - Nonlinear  Springs 


(.1)  Fixed-Base  Modes 


For  substructure  O : 


kQ  = = k2  = k ; k = k = -k2  = -k 


m(5  = - 2m  ; = m2  = m 

c ~ ? k2  k S “1 

Tc  = 1 ; <t>^  = 1 ; m0=m;  ^ * -*■  = * ] 

fa  fa  2 fa  m2  m g— * — 

Equation  (29)  yields:  do  = 2^  = - 6 ; dQ2  = do  + $ {(jlf0 

Thus  Eqs.  (15),  (30),  (32)  and  (31)  become: 


, TT n El  2 . nTT  — _ , v 

V~4_-  qSn  + Sln  V F = QSn  (1) 


NONLINEAR 

SPRING 


F = 2md(J  - k(da2  - dal) 


m da2  + k(da2  " dal>  = 0 


(iii) 


where  F = F(6,  6)  and 


« - as  - d0  - I qs„  sin  f 


sl"  2 - d0l 


(ii)  Free-Free  Component  Modes 


ka  = ka  “ k2  * k ; k = k * ~k2  = -k  ; 


ma  = ml  = 2m  ’ ma  = m2  = m ’ 

11  3 0 

*o  = ; ^a  “ 

1-2 

0 6 


0 0 


0 3/2 


tv.  = 1 


FIXED  - BASE  MODEL 


d <r=  da-| 


NONLINEAR 
fLl/-  SPRING 
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Thus  from  Eqs . (15)  and  (36): 


, tt  n El  , 2 . nir  T,  _ n 

%n  + • PA  ^Sn  + pA£  Sin  2 F QSn 

(iv) 

3mq  - F = 0 

Mal 

(v) 

6"  v + <6">l  ; v - f ' 0 

(vi) 

From  Eq . 

(33).  dQl  - qol  + ^2  > do2  - Iq!  2%2 

Thus  (v) 

- (vi)  and  2(v)  + (vi)  yield 

md"a2  - k(aal  ~ Sa2>  = 0 

(v) 

2mdal  + k(dal  ~ V - F " 0 

(vi) 

where  F = F(6,  6)  and  6 = d - d = £q 

b U 1 bn 

. nir 
sin 

--  19 


REFERENCES 


[1]  W.C.  Hurty,  "Dynamic  Analysis  of  Structural  Systems  Using  Component 
Modes",  AIAA  J.,  Vol . 3,  No.  4,  April  1965,  pp . 678-685. 

[2]  R.R.  Craig  and  M.C.C.  Bampton,  "Coupling  of  Substructures  for 
Dynamic  Analysis",  AIAA  J.,  Vol.  6,  No.  7,  July  1968,  pp.  1313-1319. 

[3]  W.A.  Benfield  and  R.F.  Hruda,  "Vibration  Analysis  of  Structures  by 
Component  Mode  Substitution",  Vol.  9,  No.  7,  July  1971,  pp.  1255-1261. 

[4]  D.  Ranlet,  F.L.  DiMaggio,  H.H.  Bleich  and  M.L.  Baron,  "Elastic 
Response  of  Submerged  Shells  with  Internally  Attached  Structures 
to  Shock  Loading",  Computers  and  Structures,  Vol.  7,  No.  3,  June 
1977,  pp.  355-364. 

[5]  F.L.  DiMaggio,  "Dynamic  Response  of  a Submerged  Elastic  Structure 
with  Elastic  Structures  Attached  to  It  by  Inelastic  Springs", 

Contract  No.  N00014-72-C-Q119 , Office  of  Naval  Research,  Weidlinger 
Associates,  Technical  Note,  August  1977. 


20  — 


APPENDIX  A - ALTERNATE:  DERIVATION  OF  DYNAMIC  EQUATIONS 


In  this  section  we  consider  the  equations  of  motion  for  the  system  as 

*) 

given  by  Benfield  and  Hruda  [3]  >except  for  the  addition  of  generalized  forces: 


M 


+ K 


vnb  \ / 'nb  l / "nb1 

When  using  fixed-base  modes  for  component  b and  free-free  modes  for 
component  a.  the  system  modal  mass  and  stiffness  matrices  are 


c cT 

M = 


m 

n 

0 " 

k 

o " 

u 

Tc  . Kc  = tcT 

5 5 

3 

0 

v 

0 

k, 

o_ 

(A-l ) 


(A-2) 


and  the  system  coordinate  transformation  matrix  is 


T5  " 


T V* 

cb  a 


0 

0 

*nb 


(A-3) 


The  generalized  forces  are 


T'  -T-  -T  T A 
()  = <t>  f + 9 f,  + 4-  T t f, 

a a a a b a cb  b 


(A-4) 


' T A 

Sib  = ^nb  fb 


where  f = known  forces  acting  on  main  structure  (S) 

3 


(A 


= known  forces  acting  on  substructure  (a) 


For  convenience  to  those  who  might  wish  to  refer  to  the  formulation 
given  in  Ref.  [3],  we  use  the  nomenclature  that  is  to  be  found  there. 

The  relationship  of  the  nomenclature  used  in  the  body  of  this  report 
to  that  used  in  this  section  is  shown  in  Table  A-l. 

*)  For  simplicity,  the  superscript  c which  refers  to  the  fixed-base  modes 
in  [3]  has  been  replaced  by  the  subscript  nj  e.g.,  fc  + f 

’b  nb 
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TABLE  A-l  NOMENCLATURE 

WA  BH, Ref.  E 3 1 *) 


Constrained  Base  Mode  Matrix 

Tc 

Tcb 

Displacements 

V dO 

Vqb 

Frequenc ies 

V“fo’“o 

Si 

3 

s 

c 

\ 

3 

Generalized  Coordinates 

Vq-f  o'% 

Vr'nb,r'b 

Generalized  Forces 

lVQfo’Qo 

VQnb’Qb 

Main  Structure 

s 

a 

Mass  and  Stiffness  Matrices 

VKo 

vkb 

Modal  Mass  Matrices 

V^fo’^o 

Ma,Mnb,Mb 

Modal  Matrices 

V»f  0^0 

VV^b 

Substructure 

0 

b 

4 


The  equations  of  the  system  as  given  by  Eq.  (A-l)  have,  with  vanishing  right- 
hand  side,  been  successfully  used  to  obtain  the  system  modes  and  frequencies,  in 
Ref.  [3].  It  should  be  noted,  however,  that,  for  complex  systems  having  many  dof, 
the  system  modal  mass  and  stiffness  matrices  would  require  the  use  of  an 
extensive  part  of  the  core  of  a digital  computer.  For  transient  or  shock 
problems  this  would  result  in  long  running  times,  so  that  a direct  use 
of  Eqs.  (A-l)  is  inadvisable.  However,  in  what  follows  we  will  show  that 
these  equations  can  be  simplified  so  that  they  do  indeed  yiel.  expressions 
which  are  identical  to  those  given  in  the  body  of  this  report. 


Defining  the  following  matrices 


allows  us  to  express  Eq.  (A-l)  in  the  following  form 


(A-5) 


*)  Symbols  not  appearing  in  [3]  are  chosen  to  be  consistent  with  the 
notation  found  there. 
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Va  + Ma“k  + ft 'C"bTcJ>  Va  + \\*bVnb^b 


+ ^(/Al  *)  <J>  5 + 4>TT  **k  4>  C . = Q 

ra  cb  b cb  a a a cb  b nb  nb  a 


M '£  . + M ,o )2r  + (f)1,  m,  T i + t/,  k,  T *4>  £ = Q , 

nb  nb  nb  nb  nb  nb  b cb  a a nb  d cb  a a nb 


As  a consequence  of  the  constraint  base  mode  definition  in  Ref.  [3], 


WA  ■ <Eb  + VclA^a  ' ° 


except  at  the  interface  coordinates.  However,  since  <j>  , = 0 at  these 

nb 


coordinates. 


^nbS/cb^a  = 0 


Also,  assuming  a lumped  mass  approximation,  we  have 


♦nbVcb  = ft,Vcb 


Upon  substituting  the  identities  of  Eqs.  (A-9)  and  (A-10)  into  Eq.  (A-7), 
we  obtain  the  final  form  of  the  equations  of  motion  for  the  fixed-base 
(constraint}  modes  of  component  b 


M £ , + M to  ^ m T £ = 0 , 

nb  nb  nb  nb  nb  nb  b cb  a a up 


Let  us  now  consider  all  but  the  first  two  terms  on  the  left-hand 
side  of  Eq.  (A-6)  and  designate  them  by  the  symbol  LHS.  Appealing  to 
the  following  identities 

*T  * 


TcbVcb  = \ + TebV 


cb 


TckU  “ TcPb*nb 

Tcb\Tcb  " Rb  + STcb  + ftb*b  + TcbVcb 

ChU  * Vnb  + TJVnb 


(A-6) 


(A-7) 


(A-8) 


(A-9) 


(A-10) 


(A-ll) 


(A-12) 
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and  rearranging  terms,  we  write 

LHS  = + (k,  + k T )$  £ + TT,  (k,  + k T ,)$  £ 

a a b bcbaa  ebb  bebaa 


(A-13) 


+ T,<rn<J)  £ + k,  <J)  , + m T 

b nb  nb  cb]  b nb  nb  b nb  nb  b cb  a s 


Now  the  equilibrium  equations  (A-ll)  can  be  written  as 


^nbj^^nb^nb  + ^b^nb^nb  + ^^cb^i^a  / ^nb 


or 


mb^nb^nb  + Wnb  + VcbV’a  = fb 

The  identities  of  Eqs.  (A-8)  and  (A-15)  allow  us  to  rewrite  Eq . (A-13)  as: 


(A-14) 


(A-15) 


LHS  = 4>T 
a 


Va5a  + (kb  + VcW’a  + VnbCnb  + Tcbfb 

Finally,  substituting  Eq.  (A-16)  into  Eq.  (A-6) , we  obtain 

l=  Ql 


(A-16) 


2 -T 

~ - 

M £ + M io  F + <p 

a a a a a a 

m.  <J>  £ + (kL  + IcT  £ + K<P  u£  u 

baa  b b cb  a a brpb  nb 

AT  A —T  — 

= <\>  f + d>f 
a a va  b 


Equations  (A-ll)  and  (A-17) , when  combined  with  appropriate  initial 
conditions,  are  a complete  mathematical  statement  of  the  problem.  They 
are  identical  to  Eqs.  (15),  (16)  and  (20),  except  for  the  fact  that  in 
the  body  of  the  report  the  forces  acting  on  the  substructure  were 
assumed  to  be  zero. 


(A-17) 


Free-Free  Modes  of  Substructure.  It  can  readily  be  shown  that  the 
equations  of  motion  for  the  system  which  appear  in  Ref.  [3],  when 
free-free  modes  of  the  substructure  are  used,  can  be  written  as  follows: 


- 


•-  'r  --u. 
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(Ma  + $IvVCa  + V^a  + Wa  + = % 


^T/v  /\  ••  aT^  — ^T/v  /\ 

Vt>  + Vb*aCa  + WbS>  = Qb 


where  the  generalized  forces 


Q = <J>  f + <p  f, 

3.  3 3 3 0 


aTa 

Qb  = Vb 


have  been  introduced.  Noting  that 


and 


rT'  \ + *KK  ‘ "b 


♦bVb  ' »£(SA + kb»b)  + *b(i:b*b  + W 


4>  € = 

3 3 D D 


Mb, 


we  rewrite  Eqs.  (A-18)  and  (A-19)  as: 

"b*beb  + W,  + wJ  ' < 


M 5 + M w2C  + <PT 
a a a a a a 


Vb  + Vb?b  - }b 


•] 

VbS>  + Wb  + wJ  ■ Qb 


(A-18) 

(A-19) 

(A-20) 


(A-21) 

(A-22) 

(A-23) 

(A-24) 


f 

I 


^ . 
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NONLINEAR  MOUNTINGS 

When  linear  elastic  system  components  are  attached  to  one  another 
by  nonlinear  mountings,  it  is  necessary  to  isolate  those  mountings  from 
all  of  the  components  of  the  system.  Thus,  only  the  forces  transmitted 
by  them,  which  are,  in  general,  functions  of  the  relative  displacements 
and  velocities,  enter  into  the  analysis. 

Fixed-Base  (Constraint)  Modes.  The  component  displacements  are 
expressed  as 


As  a consequence  of  Eq.  (A-25) , the  potential  energy  of  components 


(A-28) 
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Similarly,  we  can  express  the  kinetic  energy  as 


(A-29) 


For  simplicity,  let  us  assume  that  components  a and  b are  attached  by 
mountings  whose  properties  are  characterized  by  nonlinear  elastic  and 
linear  viscous  elements.  We  express  the  potential  energy  of  the  mountings 
as 

(PE)m  = Vm(6)  (A-30) 


and  the  viscous  dissipation  function  as 


»M  ' V6> 


(A-31) 


Using  the  generalized  form  of  Lagrange's  equations,  and  noting  that  the 

generalized  coordinates  are  £ , £ , and  6,  we  readily  obtain  the  following 

a n d 

system  of  equations 


M C + M 0)2C+$T 
a a a a ^ a 


i5a  " 6)  + (Rb  + Vcb^a  "6)  + Vnb* 


nb 


= Q (A-32) 
a 


M.C,  + M .0)^  . + 4>  T . ($  C - 6)  = Q , 
nb  nb  nb  b nb  nb  d cb  a a nb 


(A-33) 


- 6>  + <Cb  + - 5)  + Vnb5nb 

where  Q^,  are  defined  in  Eq.  (A-4) 
and 


3V, 


H 3^1 


96  + M 'Q5 


(A- 34) 


Q6  “ -fb  “ Tcbfb 


(A-35) 


Note  that  Eq.  (A-34)  is  an  expression  for  the  forces  transmitted  by  the 
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Substituting  Eq. 

M i + 

a a 

V-b  + 

where 


(A-40)  into  Eqs.  (A-38,  39)  yields: 


M w2?,  + <£TF(6,  6)  = ^Tq  + o'  = 6Tf 
a a a Ya  a o a Ya  a 

Vb«b  ~ *bF<6’  + 


(A-42) 

(A-43) 


6 = +.«.  - V’b 


(A-44) 
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